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We calculate the medium modifications of the pirir vertex in a relativistic hadronic 
framework incorporating nucleons and A(1232) isobars, and find a substantial increase 
On . of the piTTT coupling, dominated by processes where the A is excited. The coupling 

I depends significantly on the virtuality of the p, which is related to analytic properties 

^ ' of the vertex function. We analyze the general case of a non-zero three-momentum 

^ ■ of the p with respect to the nuclear medium, and evaluate the resulting widths and 

^ [ spectral strength in the transverse and longitudinal channels for the p — > tttt decay. 

' These widths are used to obtain the dilepton yields from p decays in relativistic heavy- 

■ ion collisions. 
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^ . 1 Introduction 

In recent years a lot of efforts have been undertaken in order to understand 
the properties of hadrons in nuclear medium. This challenging theoretical issue 
gains a lot of importance at the beginning of the operation of RHIC, where 
a proper and accurate inclusion of hadrons is necessary to describe the evolu- 
tion of the "hadronic soup" formed in the collision. It is commonly accepted 
that hadrons are significantly modified by nuclear matter. The arguments range 
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from simple scaling of masses |T|,@], through numerous hadronic model calcula- 
tions ||,|,|g|g[l|,|n],|l|,|l|,|3,|lilI10' QCD sum-rule techniques , 

|,0j2|j23J2§,tomodel- 



,23^ 



approaches motivated by the chiral symmetry |£] 
independent predictions based on low-density expansion and dispersion relations 
[^,^,^,0. In these calculations masses of hadrons, or widths, or both, are 
significantly changed by the presence of the medium (for a recent review see 
[^,0,^). In addition, medium may induce meson mixing absent in the vacuum 
[^,^,0,^. An indirect indication for the modification of vector mesons has 
been provided by the dilepton production measurements in relativistic heavy-ion 
collisions (CERES ||4y], HELIOS [^). The dilepton productions has been theo- 



retically studied with medium- modified vector mesons Il42|j43| , |44| , |45| , ^ , |47| , which 
helps to explain the low-mass enhancement of the dilepton yields. 



There is a lot of studies of the meson two-point functions in the literature, but only 
a few devoted to meson three-point functions. In their study of the p-meson in- 
medium spectral function Herrmann, Friman, and Norenberg have computed 
the pTTTT vertex for the p at rest in nuclear matter. They have applied a hadronic 
model with the A isobar and with non-relativistic couplings. Temperature effects 
on the pTTTT interaction have been considered by Song and Koch [J^. Urban, 
Buballa, Rapp and Wambach have extended the calculation of Ref. to 

non-zero three-momenta of the p meson and non-zero temperature. The authors of 
the present paper have analyzed the uo — > vrvr decay in nuclear medium 



Krippa has theoretically studied the effects of density on chiral mixing of meson 
three-point functions . In Ref. |Q the pTrvr coupling has been studied in quark 
matter. Otherwise, the topic of medium effects on hadronic couplings is very much 
terra incognita. Certainly, if the two-point functions are significantly altered by 
the medium, one expects that the three-point functions should also change. The 
issue is important for modeling the hadronic evolution in relativistic heavy-ion 
collisions, since the change of the hadron couplings results in altered transition 
rates between hadrons. 



Out of many meson couplings, the pTrvr vertex is especially important, since the 
p plays an essential role in hadron dynamics. The vacuum value of the coupling 
constant is large, (^p^^r — 6. Since the couplings of p and pions to nucleons and 
A isobars are also large, we expect significant medium modifications of Qpt^t^. 
As our calculation shows, this is indeed the case, with the coupling significantly 
altered already at the nuclear saturation density. The effect is mainly due to the 
interactions with the A. Our method is similar to Refs. |l7|,[TT|, with the following 
differences: We use a fully relativistic framework, with relativistic interactions. 
The leading- density approximation is used, which makes the calculation very 
simple, as no integration over the nucleon momenta is necessary. We analyze 
analyticity in the virtual p mass, which is non-trivial in our problem due to 
low-mass thresholds. These thresholds largely influence the results. In addition. 
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keeping non-zero three-momentum of the p, as in Ref. |TT[, allows us to look 
separately on the longitudinal and transverse polarizations. In our comparison 
to the CERES dilepton data we use the fire-cylinder expansion model of Refs. 
3^,p5[] and take into account the experimental cuts, which are very important 



for the detailed numerical analysis. 



The outline and the main results of the paper are as follows: In Sect. ^ we in- 
troduce our framework, i.e. the relativistic model of mesons interacting with 
nucleons and A isobars, the latter treated as Rarita-Schwinger fields. We intro- 
duce the necessary hadronic vertices and set the coupling constants. As already 
mentioned, we work at zero temperature and to the leading-order in baryon den- 
sity. In Sect. we present our results for the pmr vertex for the p meson at rest 
with respect to the medium, and point out large medium effects even at the nu- 
clear saturation density. We investigate analyticity of the coupling in the virtual 
p meson mass and point out its significance for the quantitative results. We look 
at effects of the energy-dependent width of the A and find they are not very large. 
In Sect. ^ we analyze the general case of finite three-momentum of the p with 
respect to the medium, and calculate the widths and spectral functions in the 
transverse and longitudinal channels. In Sect. ^ we apply the obtained spectral 
functions to compute the dilepton production rate from p decays for the CERES 
Pb + Au experiment. The calculation carefully includes the CERES experimental 
kinematic cuts, which is crucial for the results. We also include the effects of the 
expansion of the fire cylinder, which enter at the level of a few percent. It is 
found that spreading of the width helps to explain the low-mass enhancement of 
the dilepton spectra, but the overall normalization fails significantly short of the 
experimental data. Section ^ presents our conclusions and a discussion of several 
additional points. 



2 The framework 



Throughout this paper we choose conventionally q as the incoming momentum 
of the p carrying isospin index b and polarization vector e^, p as the outgoing 
momentum of the pion carrying isospin index a, and q — p as the outgoing mo- 
mentum of the pion carrying isospin index c. With this convention, the vacuum 
value of the pnir vertex (Feynman rule) is 

-zV,,^^.^. = g,e^''\2p^^-q>^). (1) 

The medium modifications of the coupling, calculated below, will be compared 
to Eq. (g). 
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Our calculation of the in-medium pvrvr vertex is made in the framework of a fully 
relativistic hadronic theory, where mesons interact with the nucleons and A iso- 
bars. Other in-medium calculations indicate that the leading-density approxima- 
tion is sufficient up to densities of the order of the nuclear saturation density.^ 



Also, in Ref. [|I^ the finite-temperature effects in the pinx coupling have been 
found to be moderate for temperatures up to about ~ ISOMeV, typical for rela- 
tivistic heavy-ion collisions (cf. Fig. 7(a,b) of [|T^]). This justifies the use of zero 
temperature to obtain the vertices of Fig. 1, at least as a first approximation. 

To the leading-density order only the diagrams shown in Fig. 1 contribute to the 
p ^ TTTT process. The solid lines denote i times the in-medium nucleon propagator, 
which can be decomposed in the usual way into the free and density parts : 



iG{k)=iGF{k)+iGD{k)= (2) 
= i(|^ + itln) 



1 77r 

—5{ko-Ek)e{kF-\k\) 



k"^ — mj^ + ie ' 

where k denotes the nucleon four- momentum, is the nucleon mass, Ek 



my + k^, and kp is the Fermi momentum of nuclear matter. The double line 
in the diagrams of Fig. 1 denotes i times the relativistic A propagator. 



3 3mA 3mA 



- {niA - I Fa 



(3) 

This formula corresponds to the usual Rarita-Schwinger definition ||56|j57[| , with 
the parameter choice A = —1. We have modified the denominator in Eq. (j^) in 
order to account for the finite width of the A resonance. 

Since we are interested in density effects, one of the nucleon lines in each of 
the diagrams of Fig. 1 must involve the nucleon density propagator. Go- For 
kinematic reasons, diagrams with more than one Go vanish. The wavy line in 
Fig. 1 denotes the p meson, and the dashed lines correspond to pions. All external 
particles are on mass shell. 

The meson-nucleon vertices (Feynman rules) needed for our calculation have the 
standard form 



^ An extension beyond the leading-density approximation would require accounting 
for the Fermi-sea motion, which technically leads to keeping the integration over k in 
our expressions, and, more importantly, the inclusion of correlations. In many similar 
calculations the Fermi-motion effects have been found to influence the results weakly 
at moderate densities. The inclusion of nucleon correlations is a difficult problem, ex- 
tending far beyond the present work. 
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(a) 




(b) 





(d) 





(f) 




(g) 




Fig. 1. Diagrams included in our calculation (crossed diagrams not displayed). Wavy 
lines denote the p meson, dashed lines the pions, solid lines the in-medium nucleon, 
and double lines the A. 



NN 



9a 



2mN 

; 9p9A ^fj,^ ^abc^ 

2F^ 



(4) 

(5) 
(6) 



where p is the outgoing four-momentum of the pion, q is the incoming momen- 
tum of the p meson, and a and b are the isospin indices of the pion and the p, 
respectively. We have chosen the pseudovector pion-nucleon coupling. The vertex 
(PD follows from the minimal substitution in Eq. (HI). 
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Whereas the couphngs of mesons to nucleons hsted above are well established 
and the corresponding parameters are well known (except for the long-lasting 
controversy with Hp), the relativistic couplings involving the A resonance are a 
topic of an on-going research and discussion |]57| , ^|59| , |60| . Several structures are 
possible for these vertices, and ambiguities related to the choice of the so-called 
off-shell parameters have not been resolved. For our pragmatic goal of estimating 
the size of the in-medium pirn vertex, we will adopt the following simple-minded 
and popular philosophy: all off-shell couplings are set to zero. 



The meson-nucleon-A vertices have the following form |57 



— iVpb^a^A 



where k is the nucleon four momentum, and 



(7) 

(8) 
(9) 



r^(iy) = -{W+ 2)7V, 



W = Z,Y, X. 



(10) 



The combination —{W + |) is called the off-shell parameter. According to the 
prescription stated above, we set Z = F = X = — i, such that all off-shell 
parameters vanish, and we simply have B'^^CW) = g"^. Furthermore, we take 
arbitrarily g2 = 0. The matrices in Eqs. (0-^ are the standard isospin | — >■ | 
transition matrices, given in App. 0. The vertex (P) follows from the minimum 
substitution in Eq. (0). 



For the vrAA coupling there are, according to Ref. |£9[, three possible structures: 

^K"A.A, = {G^g">'^J, + G2(7 V + P°7^)75 + GgT^I^TsT^) T^- (H) 

We drop the off-shell couplings by setting ^2 = ^3 = 0. For the pAA vertex we 
use the minimal vector current coupling [B^ and the universality prescription, 
which gives 



(12) 



The constants gi and Gi are adjusted in such a way, that in the nonrelativistic 
limit we recover the couplings ^/2{f* /m^)eijkp'^ S-j^T^ and ifA/'m'T,)PjS/^T^, re- 
spectively [|6ll , |62| , where k is the spin index of the p. The comparison, with the 
explicit form of the Rarita-Schwinger spinors and the matrices Ta and S'a (see 
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App. 0) gives Gi = \f/^/m^ and gi = \/2f*/m^. Our choice of the physical 
parameters is as follows [|62l : 

gA = 1.26, = 93MeV, = 139.6MeV, 

gp = 5.26, = 6, (13) 

^;vA = 2.12, r = 2.12, /a = 0.802. 

In addition to the mentioned ambiguities in the choice of the form of the A 
couplings p7| , |58| , p9| , |60| , there are uncertainties in the values of the coupling con- 



stants. For some constants one typically uses the quark model predictions, or 
the large-A^c arguments which relate the A couplings to the nucleon couplings 
^,0,0. Any adopted scheme should fit the values of the coupling constants to 



the available data for various processes. However, for our purpose of estimating 
the size of the medium effect on the ptttt vertex these ambiguities are not essen- 
tial. In addition, the result, as we shall see shortly, is dominated by diagram (g) 
of Fig. 1, containing only the vrA^A and irpNA couplings which are well estab- 
lished. Our results are not very sensitive to the choice of other couplings. Also, 
for simplicity of the approach and from the lack of knowledge we do not include 
any form-factors in the vertices. 

We work in the rest frame of the nuclear matter, however effort is made to write 
all expressions covariantly, which turns out to be very convenient. Our calculation 
is made in the following way. First, we evaluate the diagrams of Fig. 1. The result 
for the full p ^ vrvr amplitude has the generic form Aiacb = ^fi^acb^ where the 
vertex function is 

Ae,, = ^?p(2p^-g^), 



"^'-^ = 1 72^9^^^^' + + ck^eik^ - |k|). 



(27r)3 E, 

and A, B, and C are scalar functions depending on scalar products of the four- 
vectors q, p, and k, with k^ = E^. The occupation function is made explic- 
itly Lorentz-invariant when we write |k| = {k ■ m)^ — m^, where u is the four- 
velocity of the medium. The term with k^, upon the evaluation of the integral, can 
be in general proportional to the three Lorentz vectors present in the problem, 
namely 

where Cg, Cp, and Cu are scalar functions of p^, q"^, p ■ q, p ■ u, q ■ u, and kp. 
Contracting Eq. ( [I5| ) with p^, q^, and we obtain a set of linear algebraic 
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equations for Cq, Cp, and Cu, which can be solved. However, at the leading- 
density approximation the problem becomes even simpler. We can work in the 
rest frame of the medium, where = (1, 0, 0, 0). It is obvious that the leading- 
density approximation is equivalent to setting the three-vector k to zero in the 
functions A, B, C, and Ek appearing in the integrands of Eqs. (p!^JT5| ). Then 
/ |k|) = jPb, while higher-order terms in k generate terms with higher 

exponents of the baryon density ps- Now, with k = 0, A;'' = m^, the contraction 
of Eq. (^) with q^, p^, and gives the set of equations 

^PsmNP^C = Cpp^ + CqQ ■ p + CuP°, 

^PBmNq°C = Cpp-q + Cgq^ + Cuq'>, (16) 
^PBmiyC = Cpp ■ u + Cqq -u + Cu, 

where C is C with k = 0. Since in the general case vectors p, q, and u are linearly- 
independent, the solution of Eqs. (0) is Cp = Cg = 0, Cu = ^PbttT'nC. Thus, 
only the term proportional to u^^ in Eq. ([T5|) is present in the leading-density 
approximation. In other words, we have 

<ed = \pb{Ap'' + Bq^ + Cm^n^), (17) 

where the coefficients A, B, and C are obtained from A, B, and C by simply 
setting k = 0. 

To summarize this part, we restate the necessary steps needed to obtain the 
leading- density amplitude: The traces in diagrams of Fig. 1 are evaluated, leading 
to Eq. (|n|). Then we replace by m^u^, set k = 0, and arrive at Eq. (|T^. 
Certainly, this very simple method is general for any problem involving baryon 
loops with density-dependent nucleon propagators. In our calculation we have 
used a standard Dirac algebra package The isospin traces are evaluated in 
App.0 



3 Results for p decaying at rest 

We begin the presentation of the results with the case where the p is at rest 
with respect to nuclear matter, q = 0. In this kinematics we find B = ~^A, 
and C = 0. The fact that (7 = is reflecting the equality of the transversely 
and longitudinally polarized p meson propagators at q = 0, as will be shown 
in Sect. ^. The result B = —\A is consistent with the Ward-Takahashi identity 
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Qf^Af^ = D~^{jp) — D~^{p~q), where D^, denotes the pion propagator dressed with 
nucleon and nucleon- A bubbles [0 . With q = and pions on the mass shell, we 
obviously have -D~^(p) = D~^{p — q), and, consequently, Ap ■ q + Bq^ = 0, which 
immediately gives B = —\A. Therefore, for q = the in- medium vertex function 
is proportional to 2p^ — g^: 

<ed(q = 0) = \pBA{ci = Q){2p^ - q^). (18) 

The full formula for yl(q = 0) is very lengthy, hence we present below only the 
contribution from the bubble diagrams (f,g), which are simple. In App. we list 
the contributions to A{c{ = 0) from all diagrams in the special case of 171^^ = 0. ^ 
It turns out that the N — A bubble diagram (g) is the dominant one. We have 



A(^)(q=0) = i^,, 
9 \ 




2(771% -ml- \mNM){mN + m^ + \M) 
ni\{m% - (mA - f Fa)^ + + m^M) 



(19) 

+ (M ^ -M), 



with M denoting the mass of the p meson. For comparison, the contribution from 
the iV - AT bubble. Fig. 1(f), is 



In the limit of large m^, with ttt-a — ttt-n fixed, and Fa = 0, expression (p!9|) 
reduces to 



which agrees with non-relativistic calculations. 

In the following we shall treat M as the mass of a virtual p meson. Virtual p 
mesons are needed for the analysis of the dilepton production in Sect. ^. Analyt- 
icity of the vertex function is nontrivial in the variable M. We can see from the 
denominators of Eqs. (p!9| , |20|) , that for Fa = singularities occur at 



m A — mj^ — m 



2\ 2 



^2 ^ N . ^ (0.657GeV)^ (22) 

\ rriN J 



^ It should be stressed that there is no obvious expansion parameter in the problem. 
Expanding the amplitude in m-,^ or rriA — rriN and keeping the lowest terms does not 
lead to a good approximation. 
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0.2 0.4 0.6 0.8 1 1.2 

M [GeV] 



Fig. 2. The ratio of the effective couphng constant, Eq. (p4|), at the saturation density, 
to its vacuum value gp, plotted as a function of the virtual mass of the p, M. The solid 
and dashed lines correspond to the case with Fa = 120MeV and Fa = 0, respectively. 
The dot-dashed lines corresponds to the energy-dependent Fa of Eq. (|25|) . Arrows 
indicate positions of the singularities of Eq. ( p2[j2^ ). 

M^={ = (0.021GeV)2. (23) 

\mN ) 

Triangle diagrams of Fig. 1 also have the above singularities, and additionally 
bring in high-lying particle-antiparticle production singularities at = (2mAr)^ 
and at = (ttia + w^Af)^, which are physically not relevant. The above analytic 
structure is manifest in the numerical calculations presented below. For non-zero 
Fa the pole at ( p2D changes to a broader structure. Thus, analyticity is important 
— it immediately leads to large changes of the vertex function near the poles. 



Since for the case of q = the matter-induced vertex function (|T8D has the 
same Lorentz structure as in the vacuum, i.e. proportional to 2p^ — Q^-, it is 
convenient for our quantitative studies to introduce an effective pim coupling 
constant, defined as 



5'cfr 



gp + \pBA{q_ = 0) 



(24) 



The absolute value is taken, since with non-zero Fa the quantity y4(q = 0) is 
complex. In Fig. 2 we plot the ratio ges/gp as a function of the virtual p mass, 
M, at the saturation density, ps = Po- The dashed line corresponds to the zero- 
width A. We can clearly see the singularities of Eq. ( P^J2B| ), whose positions 
are indicated by arrows. The solid line in Fig. 2 shows the calculation with the 
vacuum value of the A width. Fa = 120MeV. In this case the pole at (^2]) changes 
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Fig. 3. Same as Fig. 2 for different values of the baryon density, ps^ with Fa depending 
on density according to parameterization (|27|). 

to a broad structure. The pole of Eq. (|^) remains, of course, unchanged. The 
considerable difference between the solid and dashed curves in the range of M 
between 0.6 and IGeV shows that the results are sensitive to the assumed value 
for Fa- We note that at low masses M, between ~ 0.07 and ~ 0.55GeV, the 
effective coupling g^.^ is lower than the vacuum value, thus the medium lowers 
the coupling, while above M ~ 0.55GeV the effect is opposite: the coupling is 
increased. Around the physical p mass, M = rup, the effective coupling is roughly 
two times larger than in the vacuum. For the width of the p —>■ tttt decay this 
means a factor of 4 enhancement, giving an in- medium width to the p of about 
600MeV at the saturation density. A similar estimate has been obtained e.g. in 
Refs. [^^. 

The dot-dashed line in Fig. 2 shows the result of a calculation with the energy- 
dependent width, parameterized as a function of the s variable as follows: 




As can be seen from Fig. 2, the effects of energy-dependent widths are not large 
in our analysis. 

It is interesting to look at the anatomy of the relative medium contribution to 
ges, i-e. of the quantity ^pBA{q = 0)/gp. ki M = nip = 776MeV and Pb = Po = 
0.17fm~'^ we find that the diagrams (a)-(g) contribute, correspondingly, —0.08, 
0.22-0.09i, -0.008 -O.OOSi, 0.10-O.lli, -0.11 + 0.63i, -0.008, and 0.62- 1.162, 
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with the total of 0.73 - 0.75i, for the case Ta = 120MeV, and -.08, 0.33, 0.001, 
0.25, -0.87, -0.008, and 2.09, with the total of 1.71, for the case Ta = 0. As 
advocated above, the largest contribution comes from the — A bubble diagram 



In Fig. 3 we display the results for different values of the baryon density, p^. In this 
study we have neglected the energy- dependence of the A width, but included the 
density dependence of Fa. The A broadens moderately at the nuclear saturation 
density ||6T| , |67[| , with Pauli blocking giving a decrease of Fa by about 40MeV, and 
absorption processes giving an increase by about 80MeV, such that the net effect 
is an increase by about 40MeV. Therefore we parameterize 

rMp.)^r.(i + 12M!XP£). pr) 

Ta Po 

Obviously, the effective coupling increases with p^, as can be seen in Fig. 3. 
A caveat is in place here. Our method can be trusted numerically only at low 
values of the baryon density, such that the leading-density approximation holds. 
On the other hand we can see from Fig. 3 that the effects are large already at the 
saturation density, and certainly the approximation breaks down at larger values 
oi Pb- Therefore our numerical results at high densities, here and in the following 
parts of the paper, have to be taken with a grain of salt and treated as indication 
of possible large effects rather than accurate numerical predictions. 

Finally, we complete our discussion of the in-medium pirn vertex at q = with 
Fig. 4, which shows the result of the calculation with fixed Fa, but with rriN and 
mA scaled down to 70% of their vacuum values. A decrease of that order at the 
saturation density is anticipated from several approaches [0^,^ . We notice that 
g^s is enhanced and shifted to lower values of M when the baryon masses are 
rescaled. 



4 Results for moving p 



When the p meson moves with a non-zero momentum q with respect to the 
medium, its propagation is different for transverse and longitudinal polarizations, 
defined by quantizing the spin along the direction of q. To analyze this effect we 
shall consider the width of the transversely and longitudinally polarized p mesons 
due to the decay into two pions. This allows us to present the result in a more 
compact form, rather than looking separately at the functions A, B, and C. The 
expression for the width for the decay p^ tt+tt", as viewed from the rest frame 
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Fig. 4. Same as Fig. 2 for the vacuum values of the mjv and (sohd hne), and for the 
values reduced to 70% (dashed line). Both curves evaluated at the baryon saturation 
density, and for Fa = 120MeV. 



of the medium, is 



r 



p — *7r7r 



j\M\\2'KyS^^\q-p~p'), (28) 



, 2qoJ (27r)32po^ i'2n)^2p'o 

where Ug is the number of spin states of the p meson, and J2s denotes the sum 
over these spin states. The division by go = V M'^ + in Eq. (^81) , rather than by 
M, accounts for the time-dilatation effect. We perform the phase-space integral 
in the rest frame of the nuclear medium, and obtain 



p — ^TTTT 



sm7- 



IP 



(b)\2 



■\M\'dj, 



(29) 



where J2b is the sum over the two possible kinematic branches, which can appear 
when the two-body decay is viewed from a frame where the decaying particle 
moves. The angle 7 is between the directions of q and p. The second branch 
appears only for |q| above a critical value. 



M JM2 - 4ml 

|q| > = Q'crit- 

2m^ 

Elementary kinematic considerations give 



(30) 



IP 



(1,2)1 



A^2|q| C0S7 ± gov - Aml{M'^ + q2 sin^ 7) 
2(M2 + q2 sin2 7) ' 
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qo 



,(1,2) 



7 



M2 + q2, 



(1,2) 

Po 



m. 



2 + (p(l,2))2^ 



d{qQ — Jml + — Jml + — 2r|q| C0S7 + q^) 



IT 



dr 

for q < gcrit 



= |p(l'2)| 



2m^|q| 



arcsm | — "^^ i_i — - | for g > gcrit 



(31) 



The transversely polarized p has two helicity states {ug = 2), with projection 
s = ±1 on the direction of q, described by polarization vectors while the 
longitudinally polarized p has one helicity state [ug = 1), with the corresponding 
projection s = 0, described by the polarization vector e'^^y An explicit calculation 

yields m,m 



^(+)^(+) 



£^(_)fc(_) —9 — u u 



{(f — q ■ u u^){q^ — q ■ u u^) 



^(0)^(0) 



Ml V 

q-q 



q ■ q — {q ■ uY 
{q^ ~ q ■ u u^){q'^ — q ■ u u'^) 
q ■ q — {q ■ uY 



T 



(32) 



Note that by summing over all polarizations one recovers the usual formula, i.e. 



q^q" 
q-q' 



(33) 



The tensors T^'^ and L'^'^ are defined with such signs as to form projection op- 
erators, I.e. , T^'^T^" = T'^", L^^L^ = L'^", and T^"'L■^ = 0. Furthermore, we 
have T^^qy = and L^'^qy = 0, which reflects current conservation, as well as 
T^'^Uy = 0. Through the use of these relation and Eqs. (|T^,^) we find that 



i-^Tp = E<;)^;4)^^ = - 

s=± 

l-^^l' = E4H4)^^ = 



2(^p + -PbA) 



(34) 
(35) 



s=0 



1 - - 1 - 

- {2{gp + -PbAYp^ + C*mNU^)L^'^{2{gp + -pBA)py + Cm^w^). 



Note that the value of the coefficient B in Eq. (0) is irrelevant for the widths. 



Equations (|3^ , |35D are used in Eq. (pQl) . 



Our numerical results are shown in Fig. 5. For simplicity, we have used here the 
constant Fa = 120MeV. We notice considerable dependence on q, as well as 
different behavior for the transverse and longitudinal cases. The transverse width 
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r^/r, 




Fig. 5. The ratio of the width for the p — > vrvr decay at the nuclear saturation density 
to its vacuum value, Fq, plotted as a function of the virtual mass of the p, M, and the 
magnitute of its three-momentum with respect to nuclear matter, |q|. Top: transverse 
polarization, bottom: longitudinal polarization. 



decreases with |q|, while the longitudinal does not. At lower values of M and |q| 
around O.SGeV the longitudinal width develops a hill, absent in the transverse 
case. 
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The quantity which enters the formula for the dilepton production (see Sect. |^) 
is the spectral function of the transverse and longitudinal p mesons, defined as[3] 



7r(M2 -m2)2 + (M2 + q2)r2,^ 



= ZTTT^^W^T^b^^^r:^^ P = T.h (36) 



where nip is the position of the pole. The results for At and Al, with nip = 
776MeV, are plotted in Fig. 6. We note that while at q = we obviously have 
At = Al, at larger values of q and at M around mp the transverse spectral 
strength becomes dominant. At first glance this may seem surprising, since in Fig. 
5 we have seen that at higher q we have much larger Ft than Fj-. However, the op- 
timum value of Fp, at which Ap has a maximum, is Fp = (M^ — m^) / \/M^~+~^. 
Lower, as well as higher values of Fp lead to a decrease in ^p. This explains 
the behavior of Fig. 6. We note that the transverse spectral strength. At, is 
concentrated along a ridge extending far into the large-q region. Thus, a proper 
description of propagation at finite and large values of q is needed for the de- 
scription of p mesons in medium. We note that our results are in qualitative 
agreement with Ref. |jTl[ (Fig. 14), with the somewhat different behavior of At, 
which reaches larger values at higher values of |q| in our approach.0 We also find 
qualitative similarity to the results of the altogether different model of Ref. [15 
(Figs. 6,7). There is a difference at larger values of |q| for At, manifest in the 
presence of the rim in our Fig. 6. The results for Al are very similar to 



We should stress here that our construction of the p spectral function accounts 
only for the process p — > tttt and ignores all other possible contributions, such 
as e.g. from the s-channel resonances, studied in Refs. |T^JT5[| , medium modifi- 
cations of the pions, etc. Such processes should be included in a more complete 
calculation. Still, the results of Sect. ^ depend only on the shape of the spectral 
functions of Fig. 6, and not on the physics leading to their form. Since the spec- 
tral functions obtained in many other approaches are quite similar to ours, the 
results obtained below can be viewed as representative to approaches containing 
the broadening of the p. 



^ The presence of a/ + here is related to the presence of 1/qo in Eq.(p9[), i.e. to 
the fact that we are using widths viewed from the rest frame of the medium. Had we 
used widths viewed from the p rest frame, we would have 1/M in Eg. (p9|) , and factors 
of M instead of \/ A-'P + in Eq. (|36|) . Of course, in both cases the resulting spectral 
functions are equal. 

^Note a factor of l/vr difference in our definition of the spectral functions compared 
to those of Ref. Q. 
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Fig. 6. The spectral strengths in the p channel at the nuclear saturation density, cor- 
responding to the widths of Fig. 5. Top: transverse polarization, bottom: longitudinal 
polarization. 

5 Dilepton production rate 

Measurements of the low-mass dilepton spectra Q^Hl have shown significant 
excess above yields from the final-state hadron decays. In this context, the prop- 
erties of vector mesons (especially of the p) in a hadronic environment become of 
particular interest, since the Vector Meson Dominance Model is commonly used 
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to make the estimates of the dilepton yields from vector-meson decays. There 
are numerous analysis of the effect in the literature, both in hydrodynamic ap- 
proaches [|2],^00,g, and in transport theories [0J7|J7|j7|j73,|78|J7|,|g . The 
dilepton-rate formula [pl| , ^j83|j8^j85[| from p meson decays can be written in a 
manifestly Lorentz covariant way as follows: 



dN 



d^xdM^ J (27r)3E^ 



p—re^e' 



-A{M,q- u,Pb{x)) fp 



q ■ u 



(37) 



where M is the invariant mass of the lepton pair, equal to the mass of the virtual 
p meson, A = 2 At -|- vAl is the spectral function including all polarizations, x is 
a space-time point, T{x) is the local value of temperature, and Eg = \/M'^ + q^. 
The quantity Fp^e+e- in Eq- (0) width for the process p — > e+e~. 



p— >e+e~ 



(38) 



where oiqed is the fine structure constant, and mg is the mass of the electron. 
Finally, the function jp in Eq. ( P7| ) is the thermal Bose-Einstein distribution of 
the p mesons. 



exp 



q ■ u — 2Pt, 
T 



(39) 



whith denoting the pion chemical potential |8^,86], incorporated in several 
works. This quantity in some sense mimics possible deviations of the system 
from the chemical equilibrium. 



In order to describe the problem as realistically as possible, we will include the 
effects of the expansion of the medium formed in a relativistic heavy-ion collision. 
The lepton pairs are formed in a fire cylinder which moves as a whole in the lab 
system with the rapidity apc- For symmetric and central collisions ajT'c' is a half of 
the projectile rapidity in the lab. In its own center-of-mass system (CM), the fire 
cylinder undergoes a hydrodynamic expansion. In the analysis of such a situation, 
it is convenient to rewrite Eq. ( pTj) in the variables suited to both the kinematics 
of the emission process and the geometry of the experimental setup. We introduce 
M_L = \J -|- q\ , the transverse mass of the dilepton pair, y''^'', the rapidity of 
the pair measured in the lab system, u^, the transverse four- velocity of the fluid 
element producing dileptons, and a'"^, the rapidity of this fluid element in the 
lab. With these variables we have 

q-u = M^ + cosh (y''^^ - a^"^) - q± ■ u^. (40) 

The velocity of the fluid element in the lab is a relativistic superposition of the 
velocity of the fire cylinder in the lab and the hydrodynamic flow considered in 
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the CM system. Thus we have 



t'_LCosh(a) 



a"" = a + apc = ^.Tctanh. v\\ + apc, u± = ^ (41) 

1 — f ^ cosh^fa) 



We note that the velocities v\\ and v± are defined now in the CM system. They 
depend on time and space coordinates. 

Next, we analyze the kinematic constraint of the CERES experiment. The exper- 
imental acceptance cuts can be included with help of the function 

«afe „ ^ _ I d^Pi±d^P2±dyidy2 (j) 6{Eg - Ep^ - Ep^)5^^\q - Pi - P2) 



$(M,|/"^^g^ 



J d'^Pi_Ld^P2±dyidy2S{Eq - Ep^ - EpJ5(3)(q - - p^) 

(42) 



where pi^2 are the momenta of the emitted electrons, yi^2 are the electron rapidi- 
ties, and is a product of step functions which enforces the experimental setup 
conditions: 2.1 = r/^m < 2/1,2 < Vmax = 2.65, > 200 MeV, and 9ee > 35mrad. 
Due to the smallness of the electron mass, we can assume here that rapidities 
and pseudorapidities of the electrons are equal. The construction of the function 
$(M, y'-"-'^, qj_) requires a numerical calculation of a two-dimensional integral of a 
function involving a product of step functions, which is very easily accomplished 
by a Monte Carlo method. With the inclusion of the experimental acceptance 
cuts, the dilepton production rate is 



dN 2M^ f d^q^ 



d^xdM At] iVmax - Vmin) J (2vr 



dy'-'^Tp^,+,-AU (43) 



One should stress here the relevance of the inclusion of the kinematic cuts for the 
obtained results. The function $ infiuences mostly the overall normalization of 
the cross section, and not so much the dependence on M. 

In order to calculate the dilepton spectrum, one has to assume a model of the 
hydrodynamic expansion of the fire cylinder. We adopt the fire-cylinder expansion 
model of Refs. It is assumed that the system is in thermal equilibrium up 



to time tmax, when freeze-out occurs, and the velocities depend on space-time in 
the following way: 

Z T 

v\\{t,z) = [vz + a^t) — , v^{t,r) = {vr + art) — , (44) 

Zmax\J') fmax\J') 



where 



,{t) = Zo + Vzt + ]^azt^, rmax{t) = rQ + Vrt+]^art^ (45) 
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are the boundaries of the system at time t. The parameters of the expansion are 



as follows 54.55 



tjnax = llfm, zo = 4.55fm, Tq = 4.6fm, (46) 
= 0.5, a, = 0.023fm-\ Vr = 0, = 0.05fm-\ (47) 

and the time dependences of the temperature and the baryon density are: 

T{t) = 210MeVexp , p^(t) = 260/V{t), (48) 

where V{t) = 2nzmax{t)r'^axi'^) is the volume of the fire cylinder at time t. For 
the time dependence of the pion chemical potential, Pn(t), we assume a linear 
rise from 20MeV at t = to 80MeV at t = t^ax 11 • 



Finally, the yield of leptons produced during the expansion is 

— — = f dt [ 27rrdr [ dz i — — — — | , (49) 
dMAri J J J [d'^x dMAri / ' ^ ^ 

where dN/ {d'^x dMArj) is given by Eq. (|43|) with all the required substitutions. 

In addition to the yields of Eq. (^9[) one usually adds the contribution from vector 
mesons which remain after freeze-out. This contribution is equal to 

f^maxytmax ) Zmaxytmax J / , \ 

dMA^ = fiM) J 2-^^ J "^'[ d^xdMA, )^ 

^max {tmax ) 

where r(M) is the full width of the p meson with virtual mass M, given by the 
formula 

with ^fpTTTT = 5.98 giving T{mp) = 150MeV. The physical interpretation of for- 
mula ( pOD is that all p mesons that remain after freeze-out, decay with the yield 
proportional to the number of mesons and the branching ratio to the dilepton 
channel. Finally, the full contribution is 

dN ^ dNi dN2 
dMAr] dMAr] dMAr]' ^ ' 

Our numerical results are shown in Fig. 7. The solid line shows the yield from 
the p decays, Eq. (0), with the vacuum spectral strength of the p, while the 



20 



o 
o 



T3 



□ 95 data 
■ 96 data 




0.7 0.8 

M [GeV] 

Fig. 7. Dilepton yields for the 158GeV/A Pb + Au CERES experiment from the p 
decays. The sohd hne is the result of the calculation with Eq. (^) and the vacuum 
p spectral function. The dashed line is obtained with the medium-modified spectral 
strength, Eq. p^). Contributions from other processes or the "cocktail" background 
are not included in the theoretical curves. 

dashed curve shows the calculation with our medium-modified spectral strength 
of Eq. (|36D . We can see that the medium effects redistribute the dilepton yields 
from higher to lower values of M. This is a typical effect of broadening of the 
p, found in previous investigations. Thus the tendency needed to explain the 
experimental data is correct. We can see by comparing the dashed line to the 
data that the calculated yield falls an order of magnitude below the data. Note 
that other processes, not included in our calculation, contribute in the explored 
region of M, for instance the Dalitz decays of mesons, or the u decays. Also, note 
that we are not including the "cocktail" contribution of decays of hadrons in our 
comparison, which would not be consistent. We have found that 15-25% of the 
model yields comes from decays after freeze-out. The effects of the expansion of 
the fire cylinder enhance the yields by a few percent. We stress that the overall 
normalization of the calculated curve is sensitive to the time-integrated volume 
of the fire cylinder, and to the hydrodynamic expansion parameters p7|. 



6 Summary and discussion 



Here are the main conclusions of our investigation: the medium effects on the pirir 
coupling are large, and dominantly come from the process where the A is excited 
in the intermediate state. At the nuclear saturation density and at physical mass 
of the p the value of the coupling is roughly doubled compared to the vacuum 
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value. The increased coupling leads directly to large widths of the p meson in 
medium. We have analyzed the resulting spectral functions for the transverse 
and longitudinal polarizations, with the result that at higher values of the three- 
momentum with respect to the medium, the transverse spectral function is much 
larger from the longitudinal one. Finally, we have applied our model to evaluate 
the dilepton production from the p decays in relativistic heavy-ion collisions. 
We confirm the well-known finding that a larger p width helps to understand 
the experimental data by redistributing strength from higher to lower invariant 
dilepton masses. However, the overall norm of the dilepton yields from the p 
decays is almost an order of magnitude too small to explain the experimental data. 
We remark that the results for the dilepton yields from p decays are only sensitive 
to the actual shape of the spectral functions and to the hydrodynamic expansion. 
In that sense they are not sensitive to detailed modelling of the dynamics, as 
long as the resulting functions At and Al are similar. Thus our results confirm 
the statement of Ref. |8^, namely that hydrodynamic models have problems in 



explaining the dilepton data unless the hydrodynamic evolution is exceedingly 
long. 

For simplicity, we have used on-shell hadronic couplings throughout our studies. 
Since the nucleon and the A in the hadronic loop of Fig. 1 can be off-shell, 
additional coupling structures may be present. In addition, sideways form factors 
could be included for the particles off-shell. Presently, this has not been done, 
again for simplicity and from the lack of knowledge as to how to introduce and 
choose these form factors. One could also include the medium modifications of the 
meson-baryon couplings, a feature advocated, e.g., in the model of Ref. pBt 



The width of the p meson picks up contributions not only from the pion-loop dia- 
gram, included in our work, but also from other processes. In particular, one can 



include the s-channel resonances, as studied e.g. in Refs. IT^JT^. Such processes 
can and should be included in a more complete and realistic calculation. Also the 
diagrams of Fig. 1 could in principle be supplied with higher resonances, at the 
expense of having more not well known parameters. 

In addition to the effects of the Fermi sea, studied in this paper, vacuum polariza- 
tion effects may influence the pTrvr coupling. To have an estimate of these effects, 
we have done a Walecka-type calculation where in the diagrams of Fig. 1 we have 
included the free nucleon propagators only. We have found a ~10% increase of 
the coupling when the nucleon mass is reduced from 939MeV to 700MeV. We 
have applied Pauli-Villars regulator with the cut-off parameter of IGeV in order 
to truncate high momenta in the loop. Similar order of the effect is found when 



the Nambu-Jona-Lasinio model is used along the lines of Ref. [5U|, and the mass 
of the quark is scaled down as expected from the medium effects. Thus, vacuum 
polarization effects on the pTrvr coupling are estimated to be less significant than 
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the Fermi-sea effects presented in this paper. 

Concerning the methods apphed, we have shown in detail how to obtain very sim- 
ply the leading-density approximation for diagrams involving loops with density- 
dependent nucleon propagators. We have also demonstrated how to carry a 
Lorentz-covariant calculation of the transverse and longitudinal spectral func- 
tions with help of the formulas of Sect. |. These methods and formulas, although 
very straightforward, are not, to our knowledge, commonly know. They can be 
useful in studies similar to ours. 



One of us (WF) thanks Dariusz Miskowiec for a discussion of the CERES exper- 
imental cuts. 



A Rarita-Schwinger spinors 



The Rarita-Schwinger spinors are defined as 

1 13 



(A.l^ 



with 



niAiEA + ^a) 



EA = \/f + m\. 
The polarization vectors are defined as 



^0 



and 



(A.2) 
(A.3) 



(A.4) 



ma(p, s) 



IEa + mA 
2mA 



\ EA+rriA 



(A.5) 



where xi^) is the two-component spinor. The spinor 'u'^(p, sa) satisfies the con- 
ditions 7^M^(p, sa) = and p^u'^{p, sa) = 0. 
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B Isospin algebra 



The isospin | ^ | transition matrices are defined through the Clebsch-Gordan 
coefficients as foliows: (|, /slT^'ji, ^3) — {^li^fill^jl^,) , with ^ and I3 denoting 
the isospin of the nucleon and A, respectively. In Cartesian basis the explicit form 
reads 



J.1 



( 



1 

■^2 









7a 



7a 
7a 



71/ 












fB.ll 



2' 



i left to right, and the rows by 



where the columns are labeled by 

h — \,\,—\,—\, top to bottom. The following useful relation holds: 



(B.2) 



The couplings of the A to the isovector proceed via the matrix defined as 



15/3 
2 \2 



Ug/illll/a). Explicitly, we find 



-'a 



'^0^0 ^ 
^010 
10^ 

V f y 







Ti 



V2, 

2 




2 " 





-1 


" 2 



0^ 



y 



I 
I 
-| 
-I J 



(B.3) 



where the columns are labeled by /g = f'i'"!'"!' right, and the rows 

by h 



3 1 



2' 2' 2' 



^, top to bottom. The conventional factor of vl5/2 in the 
definition ensures that T\ simply measures the third component of the isospin 
of the A. The spin coupling matrices used in non-relativistic calculations, 5** and 
5'X, are defined analogously and have exactly the same values as T" and T^. 



Isospin trace factors for diagrams of Fig. 1 can be now readily obtained. They 
are equal to 2, — |, — |, |, |, 2, and | times — ie"'^*, for diagrams (a),(b),...,(g). 
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respectively. Isospin indices a, b, and c have the assignment as specified on the 
text above Eq. (|l]). 



C Amplitudes for q = 



For the case q = the amplitude corresponding to the diagram (i) of Fig. 1 can 
be written as 

where G'(j) are products of coupling constants: 



_ 9a9p ^ _ V^9Af*fnNA ^ _ 9pfl 
(a) g^2' 2F^ml ' 2ml 



NA 



\mi 

ri _ 9pf^NA n — 3V2/*/7rAfA/A ^ _ 9p9a 



9pfnNA 



Gia) = (C.2) 
ml 



The formulas for A'"(j) and Z^(j) are very long in the general case, which refiects 
the presence of many terms in the Rarita-Schwinger propagator. Numerators N(^i-j 
become manageable in the formal case Fa = 0, and m.„ = 0, where we find 



A''(b) = 16 m 



4 m% + Kp m% M' 



iV(,) = -16 



{mN + mA - M) (m^ + mA + M ) 
(2 m% + rriA — ttln ttia^ + 2 mA^ — 2 (^2 mA^ + 
mN (mN - tua) {rriN + niA^ 



1 + Kp) m% mA - 2mN mA^ + (1 + i^p) "^A^) - 

-A + Kp)m% + A{l + Kp)m%mA-\ . ^ . 

^ (2 + S/tp) m^mA^ + 2/tpmA^ j 



N^a) = -20 (4 m% mA - niN niA^ - 4 7- ^ 



tmA 



m 



'AT 



7mA^ + M^ 
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/ / 

V 



N(^) = 40 



2 m%- + 4 m% m^— 
2mNmA'-7mA' 



(mjv + mA) I I 



I^A — 2 m^v ?7T.A + 2 mA*" 



mA {iTiN + ttia) (2 + mNiTiA — 2 mA^) + ^ 



7V(^) = 16 



^4 (mjv - ^a) niA^ {rriN + mA)^+ 



(C.3) 



The expressions for D^) for the case Fa = are 



D 



(a) 



D 



D 



D 



id) 



4:771% -niN M^) (-m^^ + m% M^) , 
9 tua^ {niN + niA — M ) [ttin + tua + M) (mlf — tua^ + ^tt^ — rriN 

9mA^ (-4m^ + mArM2) (^(m^ - mA^ + m/)^ - M^^ , 

27mA^ (^(m% - tua'^ + m^'^y - m% ^ , 
£)(e) = 27mA^ {niN + mA - M) {niN + mA + M) x 
(^m^ - mA^ + m^^) - M^^ , 

D(g) = 9mA^ ((m^ -W + m/)^ -m^M^^ . (C.4) 
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